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Abstract 

The magnetocaloric effect is the thermal response of a material to an external 
magnetic field. This manuscript focuses on the physics and the properties of 
materials which are commonly used for magnetic refrigeration at cryogenic 
temperatures. After a brief overview of the magnetocaloric effect and as- 
sociated thermodynamics, typical requirements on refrigerants are discussed 
from a standpoint of cooling power density optimization. Finally, a compi- 
lation of the most important properties of several common magnetocaloric 
materials is presented. 
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1. MAGNETIC REFRIGERATION 

Magnetic refrigeration is employed to provide cooling over many orders of 
magnitude of temperature, from a few hundred Picokelvin which have been 
reached with nuclear demagnetization refrigerators [1] to several hundred 
Kelvin [2, 3]. The magnetocaloric effect, which provides the basis for these 
refrigeration techniques, is discussed by Shirron [4], For low temperature 
applications, the reader’s attention is also drawn to Ambler’s and Hudson’s 
review of magnetic cooling below 1 K [5]. 
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Adiabatic Demagnetization Refrigerators (ADRs) are typically used in 
the temperature range between a few Millikelvin and a few Kelvin. In this 
temperature range, ADRs mainly compete with 3 He/ 4 He dilution refrigera- 
tors and 3 He sorption coolers; ADRs have the advantage of more efficient and 
more accurate temperature control and are thus commonly used to cool low 
temperature detectors, both in laboratory environments and on spacecraft 
(e.g. [6, 7]). The thermodynamics of ADRs are described in [8]. 

In Active Magnetic Regenerators (AMRs), the refrigerant material and 
the regenerator material are combined into one. This concept allows for 
the construction of efficient and compact refrigerators. AMRs are being 
researched for gas liquefaction (e.g. hydrogen [9]), as well as for room tem- 
perature refrigeration [10]. The working principle and the thermodynamics 
of AMRs are described in [11], 

The basic principles underlying magnetic refrigeration (with the exception 
of nuclear demagnetization, which is discussed in [12]) will be reviewed in this 
introductory section. For more details, reference is made to textbooks such 
as [13], [14], [15] or [16]. 

1.1. The Magnetocaloric Effect 

The entropy of a magnetocaloric material is the sum of three independent 
constituents, namely the magnetic entropy S m , the entropy of the lattice Si 
and the entropy of the conduction electrons S e : 

S = S m {B,T) + S l {T) + S e {T). (1) 

The magnitude of all constituents depends on the temperature T, while 
the magnitude of S m has an additional strong dependence on the external 
magnetic field density B. Si and S e can be neglected at low temperatures, 
as S m is much larger than Si and S e combined in that regime. 

On a microscopic level, the magnetocaloric effect is caused by the inter- 
action of the external magnetic field with magnetic moments in the material. 
In most cases, the magnetic moments originate from the presence of unfilled 
electron shells in transition metal or rare earth ions in the refrigerant, which 
imparts a total electronic angular momentum J on each ion. The magnetic 
moments thus have 2J + 1 possible orientations, and - assuming that the 
thermal energy ksT is greater than the splitting between the energy levels 
of these states - all will be equally occupied and will contribute 

S m = R • ln(2J + 1) (2) 
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to the magnetic entropy, where R is the gas constant (8.314 J mol” 1 K” 1 ). 

If an external magnetic field is applied, the magnetic moments will tend 
to align with the field, increasing the energy difference between the states, 
reducing the occupation of higher energy states (those align counter to the 
field), and thus reducing magnetic entropy. Using statistical mechanics, the 
field dependence of the magnetic entropy can be derived from the partition 
function. Assuming that the magnetic moments are completely independent 
and do not interact with each other, one obtains 


% = x-(coth(x)-(2J+l)-coth(x-(2J+l)))+lnf Smh(:r .‘ , (3) 

R \ sinh(x) J 

where x depends on the temperature T, the magnetic field density B, and 
properties of the refrigerant. 

For a paramagnetic material, x can be calculated with the Bohr magneton 
Hb (9.27- 10” 24 J/T), the Boltzmann constant k B (1.38 • 10” 23 J/K) and the 
Lande g-factor g (an electronic property of the magnetic refrigerant which 
will be discussed in more detail in section 2.1). 

t-iRgB 

X ~ 2k B T ' 

Equations 3 and 4 show that the magnetic entropy of a material depends 
on the ratio B/T. As B approaches 0, equation 3 will approach equation 
2, independent of temperature. In real materials, interaction between ions 
will always induce some splitting of the energy levels, causing the entropy to 
approach zero at low temperature. Although these interactions are complex, 
a rough approximation, useful in ADR design, is to replace B with an effective 
field B e s in the equations above [13]: 


= VB 2 + b 2 , 


where b is a material-dependent parameter. 

At low enough temperature, most paramagnetic materials undergo a 
phase change to an ordered state, usually ferromagnetic or antiferromagnetic 
(see section 2.2). Entropy drops sharply below the critical temperature. Fig- 
ure 1 shows the dependence of the entropy S on the magnetic field density 
B and temperature T for Ferric Ammonium Alum (a common paramagnetic 
refrigerant). The phase change can be clearly seen by the sharp change in 
slope at 26 mK. 
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Figure 1: Entropy vs. temperature curves of Ferric Ammonium Alum (FAA) for different 
values of the external field density B. The entropy is normalized to the gas constant R. 

1.2. Thermodynamics of the Magnetocaloric Effect 

The magnetocaloric effect can also be described using formalisms from 
classical thermodynamics. The total differential of the entropy S(T, B ) of an 
isobaric system is 


dS 




( 6 ) 


In an adiabatic system, dS = 0. Using the definition of the specific heat 
Cb at constant magnetic held, 


c * = r '(i 


B 


and Maxwell’s relation for the magnetization M, 


( 7 ) 
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one obtains the following relation for the adiabatic temperature change dT/dB, 
which is often used to characterize the magnitude of the magnetocaloric ef- 
fect: 


= . ( d JT\ (9) 

dB C B \8T j b ' 

Equation 9 shows that the magnetocaloric effect is large if the temper- 
ature T is high, if the specific heat C of the refrigerant is small, and if the 
magnetization changes rapidly with temperature; the latter condition is ful- 
filled around magnetic phase transitions in the refrigerant. 

Knowledge of entropy as a function of temperature and magnetic held is 
essential for the design of magnetic coolers, and the above relations are useful 
in determining entropy from measurements of refrigerant properties. Mea- 
surements of either the specific heat C(T, B) or the magnetization M(T, B) 
allow the calculation of A Sb- It is not possible to measure the three con- 
stituents of the entropy (magnetic, lattice and electronic) separately; only 
the total entropy change A Sb is accessible by experiment. 

If the specific heat C(T, B ) is measured on a sample of a magnetocaloric 
material, the entropy can be derived with the following equation: 

dT. (10) 

B 

In some cases, it is difficult to measure the specific heat to temperatures 
low enough such that C(T',B) can be easily extrapolated to T = 0. Then, 
the entropy from equation 10 is best separated into two parts, 


S(T, B ) = 


( C(V,B) 

V v 


S(T , B) — f 
Jo 


T ° fC(T', B) 


V 


dT ' 


B 


’ T 0 




where T 0 is the lowest temperature achieved in the specific heat measure- 
ments. Usually, the specific heat data at high magnetic fields can be extrap- 
olated to T = 0, because the specific heat peak due to phase transitions is 
suppressed in high fields. For many materials, the specific heat can be extrap- 
olated with a T 3 dependence for the lattice entropy and a linear dependence 
on T for the electronic specific heat. 
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If the magnetization is measured, the change in magnetic entropy induced 
by held B , A S(T,B) = S(T,B) — S(T, 0), can be calculated by integrating 
the Maxwell relation (equation 8), which gives 

ASiT ' B) = .C{w) B / dB '- (12) 

AS can then be computed by integrating the partial differentials of the 
magnetization M with respect to the temperature T from zero to the mag- 
netic held B. The entropy change AS can be calculated for both the mag- 
netization and demagnetization process. Rewriting the Maxwell relation for 
two sufficiently close temperatures T x and T 2 yields 


A S(T,B)=f o ( 


(M{T X ) - M(T 2 ) 


T\ - T 2 


■dB' = 


B' 


Ti — T; 


2 Jo 


(Af(Ti) - M(T 2 )) • dB'. (13) 


Thus, the entropy change can be obtained from the area between two mag- 
netization curves at different temperatures. Magnetization measurements 
can be performed very accurately with SQUID magnetometers. 

Note that the magnetization technique only yields AS (T, B). Full deter- 
mination of S (T, B) still requires some specihc heat data, to determine at 
least one total entropy vs. temperature curve. 

Determination of AS by these two methods usually yield identical results 
for paramagnetic or soft magnetic materials. For materials with hrst order 
phase transitions, the results obtained with the two methods sometimes differ 
from each other, however. The latent heat associated with the hrst order 
phase transition causes a discontinuity in the specihc heat data. 


1.3. ADR Cycle (Carnot Cycle) 

Adiabatic Demagnetization Refrigerators (ADRs) will be discussed here 
as an example of refrigerators employing the magnetocaloric effect. Details of 
another popular magnetic refrigeration concept, the Active Magnetic Regen- 
erator (AMR) Refrigerator, are discussed elsewhere [17, 18]. Regenerative 
thermal cycles - such as AMRs - are mainly used at temperatures above 
approximately 15 K, where high cooling power densities cannot be achieved 
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with the Carnot cycle clue to the exceedingly large contribution of the lattice 
entropy to the total entropy (for details, see section 2.3). 

ADRs use the magnetocaloric effect to lift heat from an object at temper- 
ature T co i d to a heat sink at T hot using the Carnot cycle. The cycle is shown in 
Figure 2. This entropy-temperature diagram shows various entropy curves 
for a paramagnetic refrigerant (FAA) for different external magnetic field 
densities. 

2.00 
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Figure 2: The T-S diagram of a typical ADR cooling cycle. With an ADR, the ideal 
Carnot cycle can be approximated very closely. The entropy is normalized to the gas 
constant R. 

In a typical refrigeration cycle, the paramagnetic salt is first thermally 
connected to a heat sink with temperature Thot (3 K in the example in Figure 
2); the magnetic field is then increased (from 0 to 3 T in our example). In 
the T-S diagram, this corresponds to the magnetization (1)— >(2). The heat 
sink must absorb the energy 

Qhot = T hot • (, S 1 - S 2 ) = T hot • A /Shot- (14) 
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The thermal connection to the heat sink is then open, leaving the salt 
thermally isolated. The magnetic field is then decreased, corresponding to 
adiabatic demagnetization from 3 T to 0.1 T; this corresponds to (2)— >(3) 
in Figure 2. During this step, the refrigerant cools down to T co i d - from 
3 K to 0.1 K in our example. In a real system, it is of course impossible to 
achieve perfect thermal insulation of the refrigerant. Achieving good thermal 
insulation is one of the key design drivers for refrigerant suspension systems 

[19]. 

Once the desired temperature is achieved, the rate of decrease of the 
magnetic field is controlled to maintain constant temperature, so that the 
heat absorption by the refrigerant just balances the heat flux from the load; 
this corresponds to (3)^(4) in Figure 2. During this isothermal process, the 
refrigerant absorbs heat 


Qcoid — T cold ■ (S 4 — S 3 ) — T cold ■ AS' co id. (15) 

In our example, this isothermal demagnetization is performed until the 
external field reaches zero. The refrigerant is then disconnected from the 
cold reservoir, and the magnetic field is ramped up until the refrigerant’s 
temperature corresponds to the temperature of the heat sink, (4)— >(1’). The 
cycle can then start over with the isothermal magnetization. 

Note that once heat sink temperature and maximum field are set, the 
heat absorbed per cycle is a simple function of T co \ d . 

The cooling power P of the ADR refrigeration cycle is proportional to 
the frequency / of the cycle: 

P — f ' Q cold = / ' ^cold ‘ AS^oia- (16) 

The cycle frequency / is limited by thermal resistances between the re- 
frigerant and the heat sink, and between the refrigerant and the cold load, 
respectively. Because of these resistances, entropy is generated in flowing 
heat into and out of the refrigerant, reducing the efficiency of the cycle as 
the frequency increases. The minimization of the thermal resistance between 
the refrigerant and the cold reservoir is thus key for optimizing performance; 
this is discussed in section 4. 

The Coefficient Of Performance (COP) of a thermodynamic cycle is the 
ratio of the cooling power Q co id to the work Q hot — Qcoid which needs to be 
performed to achieve the cooling: 



COP 


(17) 


Q cold 

Qhot Qcold 

in the Carnot cycle, this is equivalent to 


COP M «„ = . (18) 

r hot -L cold 

For an ideal ADR operating between 0.05 K and 2 K, the COP amounts 
to 0.026. The thermodynamic efficiency of a real cycle can then be defined 
as the ratio of the COP of the real cycle to the COP of the ideal Carnot cycle 
operating between the same two temperatures: 


V 


COP 

COPideal 


(19) 


1-4- Important Properties of Magnetocaloric Materials 

To be useful for the construction of efficient magnetic refrigeration sys- 
tems, refrigerants need to fulfill certain requirements. Yu, for instance, has 
summarized these requirements for AMRs [20]. 

To achieve a high cooling power density, the material should exhibit a 
large magnetocaloric effect in the temperature range in which the system 
operates. A large magnetocaloric effect is achieved if 


• the total angular momentum number J and the Lande factor g are 
large; 

• the refrigerant’s ordering temperature is sufficiently close to the tem- 
perature range in which the system operates to profit from a large 
dM/dT, but sufficiently far away to avoid the spontaneous entropy 
reduction brought about by the magnetic ordering. 

In addition, to allow high operating frequencies / and, in consequence, 
large cooling powers, the magnetic material should exhibit 


• small magnetic and thermal hysteresis; 

• a large electrical resistance to keep eddy current heating low; 

• a large thermal conductivity to minimize temperature gradients. 
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Furthermore, certain practical requirements exist: The refrigerant must 
be commercially available, and its cost must be reasonable. It is desirable 
that the material is non-toxic, does not corrode, and that the manufacturing 
process which is required to bring the refrigerant into the form in which it 
will be used in the system (e.g. crystal growing or sintering) is a simple one. 
In addition, it is desirable that a simple method exists to make good thermal 
contact to the refrigerant. 

Unfortunately, some of these requirements contradict each other; high 
thermal conductivity often goes hand in hand with high electrical conductiv- 
ity, for instance. 

2. OPTIMIZATION OF COOLING POWER DENSITY 

Usually, the magnetocaloric material for a specific application is chosen 
such that the cooling power density of the entire system is maximized. The 
optimization process needs to take into account both the magnet and the 
refrigerant. In most cases, the magnet and associated shielding dominates 
the mass and volume of the refrigeration system, but the magnet mass is 
proportional to the volume of the refrigerant. Thus, the usual design goal is 
the minimization of the refrigerant size. In non-continnons, or ’’single shot”, 
ADRs, cooling capacity (cooling power times hold time) sets a lower limit on 
the amount of refrigerant. Since large magnetic field changes facilitate high 
cooling power densities, superconducting magnets are usually used. However, 
because critical current drops with increasing field, the mass of a supercon- 
ducting solenoid increases faster than linearly with held, and, thus, there is 
an optimum held strength. 

Various optimization criteria for magnetocaloric materials have been pro- 
posed; Wood and Potter [21] suggest choosing refrigerants that maximize 
the product (Tkot — P co id) ■ A^oid that can be achieved with a given ther- 
modynamic cycle. Wood and Potter call this product ” refrigerant capacity”. 
They argne that in magnetic refrigeration, (Phot - Pcoid) and cannot be 

chosen independently from each other for a given magnetic held excursion, 
and that it consequently does not make sense to maximize one or the other 
individually. Also, with equation 15, 


(Phot - Pcoid) • AS'cdd = Qcoid ■ Thot _ TcoId (20) 

2 cold 

i.e. Wood’s and Potter’s criterion maximizes the product of the heat that can 
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be absorbed at T co ^ and the fractional temperature span across which the 
heat has to be pumped. Yet another interpretation of Wood’s and Potter’s 
criterion is that it maximizes the heat Qcoid that can be absorbed per COP 
point (see equation 18). 

Here, an alternative optimization criterion is proposed, namely the max- 
imization of AS'coid — Qcoid /T cc m per COP point. In cryogenic applications, 
where T lot is often much larger than T co id, the optimization criterion can be 
simplified: 


Q 


T] 


hot 


T 


cold 


cold 


T 2 

1 cold 


= AS, 


Tv 


hot 


T 


cold 


cold 


T 


A S', 


T 


hot 


cold 


cold 


T 


cold 


( 21 ) 


2.1. Ions and ion density 

An important factor in maximizing cooling power density is the use of 
materials with large zero-held magnetic entropy, so that the application of 
a magnetic held can result in a large entropy change. Equation 2 indicates 
that the total angular momentum J should thus be as large as possible. 

Hund’s rules [22, 23] can be used to estimate the total angular momentum 
J of rare earth ions. The electrons’ spin angular momentum S (not to be 
confused with the entropy S used above) and the orbital angular momentum 
L couple weakly via spin-orbit interactions; as a consequence, the two angular 
momenta are not conserved separately, but - following Hund’s 3 rd rule - are 
conserved as the total angular momentum J = \L + S\ if the electron shell is 
more than half full, and as J = \L — S\ if the electron shell is less than half 
full. 

For example, the Ce 3+ ion has a [Xe] df 1 electron configuration; the az- 
imuthal quantum number l is 3 for the f sub-shell. In accordance with Hund’s 
1 st rule, the spin of a single electron in a subshell must be parallel, i.e. 
S = 1/2. To satisfy Hund’s 2 nd rule, the orbital angular momentum must 
be maximized, and thus L — 3 (L — 'Y^rrii, where rrq is the magnetic quan- 
tum number and ranges from — l to /). The total angular momentum of an 
isolated Ce 3+ ion is finally given by Hund’s 3 rd rule, J = \L — S\ = 5/2. 
However, interactions with the crystal held, the electric held due to neigh- 
boring ions, can lead to partial lifting of degeneracies of ground multiplcts. 
The Kramers degeneracy theorem [24] states that systems with an odd to- 
tal number of electrons remain at least doubly degenerate in the absence of 
magnetic helds. These pairs of degenerate states are called Kramers dou- 
blets. For instance, in Cerium Magnesium Nitrate (CMN), a paramagnetic 
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salt commonly used for very low temperature cooling and thermometry, the 
J = 5/2 ground multiplet of the Ce 3+ ions is split into three Kramers dou- 
blets [25]. At low temperature, only the lowest doublet is populated, and 
thus CMN is best described as a J = 1/2 material. 

Hund’s rules predict experimental results very well for rare earths, as their 
partially filled 4f shells lie deep within the ion, and have weaker interaction 
with the crystal field. For the iron group transition metals, or 3d ions, 
Hund’s rules do not yield satisfactory results, and molecular field theory 
must be used. The failure of Hund’s rules is mainly due to interactions 
between the electrons in the 3d sub-shell and the surrounding crystal field, 
which dominate spin-orbit interactions [26]. 

Experimental data suggests that 3d systems assume a ground state such 
that J = S, an effect which is called ’’orbital quenching” [26]. Examples for 
transition metals which have been used successfully in magnetic refrigeration 
include the ions Fe 3+ in the paramagnetic salt FAA (J = 5/2), Cr 3+ in CPA 
(J = 4/2), and Mn 2+ in MAS (J = 5/2). 

Gadolinium is a special case. In Gadolinium, the 4f sub-shell is half full, 
and the orbital angular momentum L vanishes as a result of Hund’s I s * and 
2 nd rules. Hence J = S = 7/2 , the largest of any ion [27]. Moreover, there 
is no distortion of the spherical 4f charge density, and consequently little 
interaction with the crystal field. Thus, Gadolinium compounds typically 
have the magnetic entropy S m = R ■ ln(8) predicted for a fully degenerate 
J = 7/2 shell, and are often used in magnetic refrigeration. 

The contribution of spin angular momenta to the magnetic moment of an 
atomic electron is characterized by the Lande g-factor g, which is defined as 


1) -L{L + 1) 

1 ) 5 

equation 22 does not hold if spin-orbit interactions are not the main type 
of coupling. Thus, in the iron group transition metal ions, g ~ 2. For a 
magnetic refrigerant, the Lande g-factor should be as large as possible to 
maximize cooling power density. 

Dysprosium compounds have very high g-factors, which can also exhibit 
spacial anisotropy. Dysprosium Gallium Garnet (DGG), for instance, has 
a g-factor of 8 and Dysprosium Aluminum Oxide (DOA) a g-factor of 14; 
for comparison, CMN, FAA, CPA and MAS all have a g-factor of 2. In 
addition, the net magnetic moment of DOA varies greatly with the crystal 
orientation; it is almost one order of magnitude larger along the [010] axis 



, J(J+l) + S(S + 
®~ 1+ 2 JpV 
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than along the [001] axis. Kuzmin [28] proposed to use this anisotropy to 
create a refrigerator in which the refrigerant simply rotates in a constant 
magnetic field. 

Naturally, materials with greater density of magnetic ions have higher 
cooling power density. The ion densities for several common refrigerants 
are listed in Tables 1 and 2. A higher concentration means a shorter dis- 
tance between the magnetic ions, and hence stronger interaction and higher 
ordering temperature. The correlation between ion density and ordering 
temperature is illustrated in Figure 3. The correlation is not perfect because 
non-magnetic ions and the type of bond influence the strength of the interac- 
tions. For example, metals generally have higher ordering temperatures than 
ionic compounds because of the Ruderman-Kittel-Kasuya-Yoshida (RKKY) 
interaction [29, 30, 31]. In a metal, localized magnetic moments spin-polarize 
the electron gas, which in turn then interacts with other magnetic moments. 
The resulting indirect interaction between magnetic moments gives rise to a 
higher ordering temperature. 

Thus, more dilute salts must be used for cooling to lower temperatures. 
Well known compounds include the hydrated salts CMN, CPA, MAS and 
FAA, which have often been used at to achieve temperatures below 100 mK. 
Gadolinium chloride elpasolite (GCE, Cs 2 NaGdCl 6 ) is an attractive candi- 
date material because of its its low ordering temperature of 35 mK, in spite of 
its rather high ion density [32], Unfortunately, Cs 2 NaGdCl 6 is deliquescent 
and consequently difficult to handle, and has not yet been used in a prac- 
tical ADR. The Gadolinium fluoride elpasolite (Cs 2 NaGdF 6 ) is hydrophobic 
and expected to have magnetic properties which are similar to the ones of 
Cs 2 NaGdCl 6 . At higher temperatures, fluorides and oxides (such as GGG, 
GLF or GdF) are often employed as refrigerants. High-spin molecular clus- 
ters such as M 1112 or Fes have recently been suggested for helium liquefaction 

[33]. 


2.2. Magnetic Ordering 

The minimum temperature which can theoretically be achieved by adi- 
abatic demagnetization of a magnetocaloric material is determined by the 
material’s ordering temperature, where a transition occurs and the magnetic 
moments align spontaneously. This results in a reduction of the magnetic 
entropy to zero over a very small temperature range, or - in other words - 
a small entropy change over a large range of magnetic field densities. The 
temperature at which this ordering occurs is called Curie temperature T(< if 
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Ion Density [ 10 21 cm 3 ] 


Figure 3: Interactions between paramagnetic ions are largely proportional to a refrigerant’s 
ion density; larger interactions between ions result in higher ordering temperatures and 
limit the minimum temperature at which a material can be used for refrigeration. 

the magnetocaloric material undergoes ordering from the paramagnetic to 
the ferro- or ferrimagnetic state, and Neel temperature T N if the ordering 
goes from the paramagnetic to the antiferromagnetic state. 

In accordance with equation 9, the magnetocaloric effect is large when 
the derivative of the magnetization with respect to temperature [dMj d T ) b 
is large; this is the case in the vicinity of magnetic order-disorder phase tran- 
sitions. In magnetocaloric materials, ferromagnetic interaction is preferable 
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to antiferromagnetic interaction, as the applied external magnetic field is 
enhanced by the internal field. 

Some paramagnetic materials, such as Gadolinium Gallium Garnet (GGG, 
see section 3.2), do not undergo a transition to an ordered phase. In GGG, no 
long range magnetic order has been observed down to 25 mK in the absence 
of external magnetic fields [34], This is attributed to geometric frustration 
of the spins [35]; the magnetic ions are positioned on well-ordered points in 
the lattice structure such that interactions between the magnetic moments 
preclude the existence of one single preferred spin orientation. Geometric 
frustration thus leads to a large number of possible ground states and sup- 
presses magnetic ordering. The temperature below which frustration induced 
spin freezing occurs in GGG is estimated to be 130 mK. In the presence of a 
magnetic field, GGG undergoes antiferromagnetic ordering (at a field of 1 T, 
the Neel temperature T,v is 380 mK) [36]. Due to the presence of a large 
number of soft modes, the magnetocaloric effect is enhanced in geometrically 
frustrated systems [37]. 

2.3. Lattice Entropy 

The lattice entropy Si can be calculated from the lattice specific heat Q, 

S,(T) = £ Odd . dT >. (23) 

Assuming that the temperature T is small compared to the Debye tem- 
perature 9d, 


rT rpt 3 

Si(T) oc jf — ■ dT' oc T 3 . (24) 

Due to the lattice entropy’s T 3 dependence, the contribution of the lattice 
to the total entropy (see equation 1) cannot be neglected at higher temper- 
atures; for typical refrigerants, this occurs at around 15 K. 

In a magnetic refrigerator, the magnet sets practical limits to the max- 
imum achievable variation of the magnetic field. A non-negligible lattice 
entropy Si reduces the contribution of the magnetic entropy S m to the total 
entropy S = S m + Si + S e ; in that case, a given change of the external mag- 
netic field decreases the total entropy less than would be the case if the total 
entropy consisted only of the magnetic term (S = S m ). As a consequence, 
the magnetic field needs to be ramped down more to achieve a certain base 
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temperature (step (2)— >(3) in the description of the cycle in section 1.3), and 
less field is available for the isothermal absorption of Qcoid (step (3)— >(4)). 
In other words, the lattice entropy adds a significant slope to the lines of 
constant field in Figure 2, which greatly reduces the temperature change 
achievable by adiabatic demagnetization from a point at a given held and 
temperature in this region. 

As a consequence, the lattice entropy should be as small as possible to 
maximize the cooling power density. At temperatures below 15 K, where 
the lattice entropy is relatively small, the Carnot cycle can be employed 
effectively and paramagnetic salts and rare-earth oxides with low magnetic 
ordering temperatures are commonly used. 

In some materials, the entropy of the lattice can even contribute to the 
cooling power. This behavior was first observed in 1997 in the compound 
GdsSi 2 Ge 2 , and is called ’’Giant Magnetocaloric Effect” [38]. The giant mag- 
netocaloric effect is attributed to field-induced structural transformations of 
the crystal. When the external held changes, the entropy difference between 
the two structures (which is a change in the lattice entropy) adds to the mag- 
netic entropy change, and amplifies the conventional magnetocaloric effect 

[39]. 

3. MATERIAL DATA 

A vast number of materials has been studied with respect to their suit- 
ability for magnetic refrigeration. While paramagnetic salts with operating 
temperatures of a few Kelvin and below were investigated intensively in the 
mid-twentieth century, the research focus has now largely shifted towards 
rare-earth compounds for refrigeration between 15 K and room temperature. 

3.1. Hydrated Paramagnetic Salts 

Hydrated paramagnetic salts contain ions with only partially filled elec- 
tron shells. The large distance between the ions (which is required for a 
low ordering temperature) is due to the large amount of water of crystal- 
lization in the salt. In the following, the paramagnetic salts which are most 
commonly used for refrigeration are described, sorted by increasing ordering 
temperature. 

Cerium Magnesium Nitrate (CMN) is a popular magnetic refrigerant be- 
cause of its low ordering temperature of just below 2 mK [40]. The crystals 
contain a large amount of water, which assures a large distance between the 
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magnetic Ce 3+ ions (about 1 nm). The ions’ distance from each other leads 
to weak interactions between them [41] , contributing to a small internal mag- 
netic field b (cf. equation 5) and a low T c . The heat capacity of CMN has 
been measured between 0.6 mK and 4.2 K by Fisher et ah [42], Crystal 
growth is complicated because the solubility of CMN is a strong function 
of temperature so care is required during growth and properties may vary 
between samples produced in different environments [43]. In experiments 
where it can be in direct contact with liquid helium, CMN is often used 
in powdered form, which affords an advantage for heat transfer because of 
the larger surface area. The powdered form is especially popular with 3 He 
experiments because CMN has an anomalously low Kapitza resistance with 
3 He [44, 45]. However, the susceptibility of CMN is strongly anisotropic 
so that individual crystallites oriented in different directions with respect to 
the magnetic field can cool to dramatically different temperatures. Establish- 
ment of thermal equilibrium through grain to grain heat transfer is inefficient 
and has the same effect as a relatively large but decreasing heat leak [43]. 
CMN is also often used for the construction of magnetic thermometers in 
which the magnetic susceptibility of a material is measured to determine the 
absolute temperature in accordance with the Curie- Weiss law. CMN (CAS 
number 13550-46-4) is currently difficult to obtain commercially. At the time 
of writing, no industrial suppliers could be found. 

Lanthanum Cerium Magnesium Nitrate (LCMN) is a modified version of 
CMN where some of the magnetic Ce 3+ ions are replaced by non-magnetic 
Lanthanum ions. This reduces the internal magnetic held b further and thus 
the magnetic ordering and base temperatures. Temperatures below 1 mK 
were achieved in demagnetization experiments of compressed LCMN powder 
samples with 90 percent of the ion sites occupied by La rather than Ce atoms. 
However, due to contributions to the internal magnetic held b aside from the 
magnetic dipole interaction of the magnetic ions, dilution cannot eliminate 
the internal magnetic held completely [46]. 

Chromium Potassium Alum (CPA) has a long history of use in magnetic 
cooling experiments with its use dating back to 1933 [47] . Its heat capacity 
was measured between 30 mK and 0.8 K by Vilches and Wheatley [48, 49] 
and between 1 K and 4.2 K by Kapadnis [50]. It has a relatively low ordering 
temperature of 10 mK [13], and its octahedral crystals can be easily grown 
to large sizes [43]. Unlike CMN and LCMN, its magnetic susceptibility is 
isotropic so that it can be used in powder form with no less cooling power 
than in its crystalline form [43]. As with many hydrated salts, CPA loses 
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its water of crystallization quite easily and so must be sealed well for use 
in vacuum. It does not attack copper, and so can be grown on a copper, 
rather than gold, thermal bus. CPA is readily availably as small crystals 
from several suppliers of laboratory chemicals at a comparatively low price 
(CAS RN 7788-99-0). 

Ferric Ammojiium Alum (FAA) is one of the most commonly used para- 
magnetic refrigerants. The specific heat of FAA was measured by Vilches 
and Wheatley [48], who also observed magnetic ordering at 26 mK. FAA is 
readily availably in crystalline form (CAS RN 7783-83-7). FAA readily dis- 
solves copper, and has even been known to dissolve away the copper wires 
inside gold plating. As a consequence, in ADR salt pills using FAA, the 
crystals are grown around bundles of gold wires, and the hermetic enclosures 
are typically constructed from 316 stainless steel. FAA breaks down when 
subjected to temperatures in excess of approximately 35 °C, so one must 
be careful to store FAA salt pills in a temperature-controlled environment, 
particularly in warm climates. 

Manganese Ammonium Sulfate (MAS) has a magnetic ordering tempera- 
ture of 170 mK [13] so is more popular cooling samples in higher temperature 
ranges. It is part of the family of Tutton salts, which contain two different 
cations similar to the alums, but with a slightly different chemical formula. 
Like CPA and FAA, its susceptibility is isotropic so it can be used in pow- 
dered form with little loss of cooling power. Its entropy has been measured 
by Vilches et al. [48, 49] down to 90 mK. 

3.2. Rare-Earth Oxides and Fluorides 

Due to their chemical stability and high magnetic ion concentration, rare- 
earth oxides have been studied intensively. Rare earth garnets are used in 
hydrogen liquefaction applications as they do not hydrogenate. Barclay [51] 
gives a comprehensive compilation of various Gd compounds. 

Gadolinium Gallium Garnet (GGG, Gd 3 Ga 5 0i2, CAS RN 12024-36-1) 
has been used as the working material in magnetic refrigerators in the tem- 
perature range from 2.0 to 20 K for the production of superfluid helium and 
for liquefying helium [52], Good quality crystals with high thermal conduc- 
tivity and chemical stability can be produced at reasonable cost. The entropy 
change of GGG becomes small at high temperatures, however. Entropy val- 
ues for GGG were published by many authors, for instance by Fisher et al. 
[53]. The g-factor of GGG is isotropic and amounts to 2. 
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Dysprosium Aluminum Garnet (DAG, Dy 3 Al 5 0i 2 ) is a magnetocaloric 
material which facilitates higher cooling power densities than GGG; this is 
due to a strong anisotropic magnetic interaction of the Dy ion. The g-factors 
of single DAG crystal are 10.6, 8.8 and 6.5 for the <111>, <110> and <100> 
directions, respectively. 

Figure 4 1 shows entropy diagrams for single crystals of GGG and DAG 
with the magnetic field in the <111> direction [54], The zero-field saturation 
entropy is higher for GGG than for DAG, as the Gd ion in the GGG has 
a magnetic moment of J = 7/2, while the degeneracy of the large room 
temperature magnetic moment of the Dy ion in DAG is lifted by the crystal 
field at low temperatures, leaving only an effective magnetic moments of 
J = 1/2 due to a Kramers doublet. As the figure shows, at temperatures 
above 10 K, DAG provides a larger entropy change than GGG because of 
its larger g-factor, which directly multiplies magnetic held in the entropy 
function (equations 3 and 4). 

When polycrystalline DAG is used as a refrigerant, the strong anisotropy 
of Dy must be taken into account. In polycrystalline DAG, the crystals 
are randomly oriented, and the anisotropic effect is averaged. The entropy 
change of polycrystalline DAG roughly decreases to the level of the single 
crystal along the <110> direction. 

The thermal conductivity of many single crystal garnets - including the 
two refrigerants discussed above, GGG and DAG - was measured by Slack 
and Oliver [55]. In general, the thermal conductivity of GGG is typical for 
single crystal dielectrics, with a strong peak in the 10 to 30 K range. The 
peak for DAG is smaller and occurs at lower temperature, possibly because 
of phonon scattering produced by low-lying electronic levels of the Dy ion. 

Gadolinium, Dysprosium Aluminum Garnet (GDAG), i.e. DAG with a 
partial substitution of Dy with Gd, has been studied to solve the anisotropy 
problem of DAG [56, 57]. Figure 5 shows the entropy change for the poly- 
crystalline (Gd iT Dy 1 _ a ,) 3 Al 5 0 i 2 samples (x = 0, 0.1 and 0.2) and monocrys- 
talline DAG in the magnetic field of 6 T [54], The substitution of Dy with 
Gd increases the magnetic entropy change and therefore the total entropy, 


1 In Figure 4 and most of the following plots, the material data is given per unit mass, 
as this value can be measured reliably. The entropy and specific heat per unit volume 
depend on the mass density of the material, which is influenced by the manufacturing 
process and may vary between samples. In the following, entropy per unit volume will 
only be shown if the relative density of the sample was well known. 
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Figure 4: Entropy diagram of single crystals of GGG and DAG for the magnetic fields of 
0, 2, 4 and 6 T [54], 


especially in the liquid hydrogen temperature range of about 20 K. 

For practical purposes, it is important that a magnetocaloric material 
is available in various shapes such as plates and small spheres. As a con- 
sequence, polycrystalline rectangular solid and spherical particles of GDAG 
were made and tested to confirm that the production process did not in- 
fluence the magnetic properties in a negative way [58, 59]. Figure 6 shows 
GDAG shaped into a 10 cm size plate and small spheres with 0.4 mm di- 
ameter, with 99.5 percent relative density and a Vickers hardness of 500. 
The diameters found in the spherical material stay within 10 percent of the 
average value of 0.4 mm. 

Gadolinium Gallium Iron Garnet (GGIG) is an Fe-modihed Gd-Ga nano- 
composite garnet. It was demonstrated that when iron is added to the garnets 
in substitution for some of the gallium atoms, the resulting magnetocaloric 
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Figure 5: Entropy change of polycrystalline (Gd x Dyi_ x )3Al 5 0i2 samples for x = 0, 0.1 
and 0.2 in comparison with monocrystalline DAG [54]. 



Figure 6: Photographs of a rectangular solid and spherical particles of GDAG. 


response is increased in comparison to the garnets without the iron addi- 
tion [60, 61, 62, 63]. The magnetic entropy changes AS of GGG and GGIG 
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(Gd 3 (Gai_ a; Fe x ) 50 i 2 ) with x = 0.125, 0.25, 0.4 and 0.5 induced by an exter- 
nal field of 5 T are shown in Figure 7 [63]. 



Temperature (K) 

Figure 7: Magnetic entropy change AS of GGG and GGIG (Gd 3 (Gai_ s Fe x ) 50 i 2 ) with 
x = 0.125, 0.25, 0.4 and 0.5 induced by an external field of 5 T. AS was calculated from 
magnetization data [63]. 

Figure 8 shows the entropy diagram for GGIG with the comparatively 
high iron content of x = 0.4 [63] for various magnetic fields. The entropy 
of GGIG with x = 0.4 increases smoothly at low temperature, and has a 
larger slope than the entropy of GGG at high temperature. This tendency 
becomes more pronounced with increasing iron content. GGIG with a high 
iron content is thus unsuitable for refrigerators using the Carnot cycle because 
the adiabatic temperature change is small. The most suitable cycles for 
GGIG with high iron content are regenerative thermal cycles such as Stirling, 
Ericsson, and the Active Magnetic Regenerator (AMR) cycle. 
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Figure 8: Entropy diagram for GGIG with x = 0.4 in 0, 2, 4, and 5 T [63]. 


Gadolinium Dysprosium Gallium Iron Garnet (GDGIG, chemical formula 
(Gdo.5Dy 0 .5)3(Gao.875Feo.i25)50i2) was studied in polycrystalline form [64], 
The substitution of Gd with Dy in GGIG has enhanced the entropy change 
in small magnetic fields and at high temperatures in a similar way as the 
same substitution in GGG. 

Dysprosium Gallium Garnet (DGG) is part of the family of antiferromag- 
netic rare earth garnets. Unlike Gd based garnets, which are only effective 
refrigerants below 15 K because they have little magnetic entropy change 
above that temperature, Dy still has a large magnetic entropy above 15 K 
and so Dy based garnets like DGG can be used up to at least 30 K [65]. DGG 
has a magnetic ordering temperature around 0.6 K, so it is most useful as a 
high temperature refrigerant [66] . The salt can be grown using the Czochral- 
ski process with the starting materials Dy 203 and Ga 203 , which are both 
readily available from commercial chemical suppliers [65]. DGG has been the 
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subject of much study. Its heat capacity was measured by Tomokiyo et al. 
and Kimura et al. [65, 67]. Its thermal conductivity has been measured by 
[68]. The thermal conductivity of DGG features an interesting dependence 
on the magnetic field, shown in Figure 9 [69]. Strong phonon scattering is 
induced by the magnetic field and changes the thermal conductivity by up 
to a factor of 70. 



Figure 9: Magnetic field dependence of the thermal conductivity of DGG [69] 

Gadolinium Lithium Fluoride (GLF, GdLiF \) is an attractive material in 
that it has an ion desity somewhat higher than that of GGG and a transition 
temperature below 0.3 K [70, 71, 72], Unfortunately, absolute values of the 
entropy as a function of temperature and field have not been published for 
this material. Figure 10 shows the magnetic entropy change of polycrystalline 
GLF, compared with that of GGG. The entropy change of GLF is 20 to 
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60 percent larger than the entropy change of GGG. A number of groups have 
successfully grown large crystals of this material, although usually doped 
with other rare earth ions for use in lasers. 



Temperature (K) 

Figure 10: Entropy change of a GLF (GdLiFA polycrystalline sample and a GGG 
(Gd 3 Ga. 50 i 2 ) single crystal sample as functions of the magnetic field and temperature 
[71]. The mass density for GLF is 5.34 g/cm 3 , and for GGG 7.08 g/cm 3 . 

Gadolinium Fluoride (GdF 3 ) has a gadolinium ion density 65.6 percent 
that of pure gadolinium, which among dielectric compounds is second only 
to gadolinium oxide. Its ordering temperature, 1.25 K, is much lower than 
that of Gd 2 0 3 (3.8 K). Full entropy-temperature curves are presented for 
fields up to 5 T in DiPirro [72], Unfortunately, GdF 3 appears to undergo 
a phase change, with significant volume change, upon cooling from melting 
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temperature, which makes the growth of large crystals very difficult. How- 
ever, there is some evidence that this phase change is clue to the presence 
of oxide impurities, and a few groups have successfully grown large single 
crystals by taking great care to eliminate these impurities [73]. 

3.3. Laves Phase Materials 

Ferromagnetic materials are required in order to achieve magnetic entropy 
changes AS which are sufficiently large for refrigeration above 15 K, where 
the thermal energy of the phonons becomes large compared to the energies 
of the Zeeman levels in an external magnetic field of several Tesla. Espe- 
cially in regenerative magnetic refrigeration, it is essential that the ordering 
temperature of the regenerator material is closely matched to the refrigera- 
tor’s operating temperature range. As a consequence, a number of different 
magnetocaloric materials is required. Several studies have been performed 
on magnetic refrigerants in the temperature range above 20 K [74], In this 
section, we focus on Laves phase materials. 

Because of their large magnetic moments of about 10 Hb P er rare earth 
ion, RM 2 (R: rare earth; M: Al, Ni, Co) Laves phase rare earth intermetallic 
compounds with a cubic MgCu 2 structure have been studied and were found 
to have large magnetocaloric potentials and suitable transition temperatures 
for hydrogen liquefaction [75, 76]. The transition temperatures can be ad- 
justed by changing the rare earth element R and/or the metal M; the nature 
of the phase transition differs among these compounds. The magnetic transi- 
tion temperatures in the RM 2 system have a strong relation to the de Gennes 
factor G of the lanthanide ion: 

G — (g — l) 2 ■ J ■ (J + 1). (25) 

In Figure 11, the magnetic phase transition temperatures are shown as 
functions of the de Gennes factor for various RNi 2 , RA1 2 , and RCo 2 com- 
pounds [75, 77]. The transition temperature T c of a single lanthanide com- 
pound scales almost linearly with the de Gennes factor. In compounds with 
two lanthanide ions, the partial substitution of one magnetic lanthanide with 
another one modifies the phase transition temperature, which is also shown 
in Figure 11. Changing the ratio of the rare earth elements can thus be used 
to tune the properties of a magnetic refrigerant to a particular temperature 
range. 

RNi 2 and RA1 2 compounds have a second order phase transition from 
the paramagnetic to the ferromagnetic state. The specific heat of ErAl 2 is 
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Figure 11: Magnetic phase transition temperatures for various RM 2 (M: Ni, Al, Co) 
compounds as a function of the de Gennes factor. 


shown in Figure 12. RC 02 compounds (R: Dy, Ho, and Er) feature intriguing 
magnetic properties [78] that result from interactions of localized 4f elec- 
trons and the itinerant 3d electrons of Co. Dy-, Ho- and Er-based RC 02 
compounds exhibit first-order paramagnetic-ferrimagnetic phase transitions, 
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while the others show second-order phase transitions from para- to ferromag- 
netic states. Above the transition temperature, a field-induced metamagnetic 
transition occurs. 



Temperature (K) 

Figure 12: The specific heat of ErAR 

It has been shown that the substitution of Dy, Ho or Er with nonmagnetic 
Y ions (which yields R 1 _ 3 .Ya.C 02 ), and the substitution of Ho with Tb (which 
yields Hoi_ x Tba;Co 2 ) can be used to control the transition temperature [79, 
80]; however the substitution results in a change of the phase transition from 
first order to second order [81]. 

The magnetic transition temperature of RM 2 compounds can also be 
tuned by the substitution of M. For example, the magnetocaloric effect of 
Ho(Coi_a;Nia ;)2 was studied by Wada et al. [82], 
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The effect of Er substitution in Er a ,Hoi_ :E Co2 and Er3.Dyx_3.C02 has been 
studied in detail as the phase transition remains of first order in these sys- 
tems. Figure 13 shows the specific heat of (Er 0 .53Hoo.47)Co2 as an example. 



Figure 13: The specific heat of (Er 0 . 53 Hoo. 47 )Co 2 [77]. 

Figure 14 shows the isothermal entropy change AS of several RM 2 com- 
pounds. AS has a peak around T c and gradually decreases both above and 
below the transition temperature. AS of these compounds increases with 
increasing field and shows no saturation at least up to 5 T. Figure 14 illus- 
trates how all RC02 compounds clearly show characteristics of a first-order 
transition. AS has a trapezoid shape formed by a sharp onset at T c and a 
shoulder-like drop with increasing temperature. The maximum value of AS 
changes only slightly when the magnetic field is increased above 3 T. The 
temperature of the peak of AS moves to higher temperatures with increas- 
ing magnetic field, since the metamagnetic transition temperatures increases 
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with increasing magnetic field. As a consequence, refrigerators using RC 02 
materials have a larger cooling power density, but a narrower range of oper- 
ating temperatures compared to RA1 2 and RNi 2 . 

Due to the complicated temperature and field dependence of the entropies 
of RM 2 compounds, numeric thermal modeling and cycle simulation are re- 
quired to design and operate efficient AMRs with these materials. 
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Figure 14: Entropy changes AS of various RM 2 materials in an external field of 5 T. AS 
was obtained from magnetization measurements. 


In magnetic refrigerators, high cycle frequencies can be achieved either 
by using refrigerants with a high thermal conductivity, or by shaping the 
refrigerant such that the thermal penetration depth is on the order of the 
characteristic dimension of the refrigerant. Figure 15 shows the thermal con- 
ductivity of cast and sintered DyAl 2 samples, which is comparable to the 
rather poor thermal conductivity of stainless steel. As a consequence, shapes 
such as small spheres or thin sheet material are typically used for the con- 
struction of regenerators. Centrifugal atomization is used for HoA 1 2 , DyAl 2 , 
and GdNi 2 . Two representative examples of spherical compounds, GdNi 2 and 
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DyAl 2 , are shown in Figure 16. Pellets of these materials exhibit sharp mag- 
netic transitions and have good potential for use in magnetic refrigeration 

[83]. 
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Figure 15: The thermal conductivity of cast and sintered DyAl 2 samples. 


4. REFRIGERANT STAGE DESIGN CONSIDERATIONS 

The choice of refrigerant has a strong impact on the details of the de- 
sign of its enclosure, and how thermal contact is made to the refrigerant. 
McCammon and Shirron discuss the intricacies of constructing salt pills in 
[84], Hydrated refrigerants, for instance, need to be sealed in a hermetic 
container to avoid de- hydration in the insulation vacuum of a cryostat. An- 
other important point is material compatibility, especially with respect to 
corrosion. 

In the following, another important design aspect - namely the thermal 
contact to the refrigerant - will be discussed in more detail, especially in 
view of its impact on a system’s cooling power density. A good thermal 
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Figure 16: Photograph of spherical RM 2 compounds [83]. Left: GdNi 2 with a particle size 
of 355-500 pm. Right: DyAl 2 with a particle size of 355-500 pm. 


contact between the refrigerant and the sample to be cooled is required to 
achieve high thermodynamic efficiency, both in steady state (the lower the 
temperature gradient between the refrigerant and the sample, the higher the 
COP) and during transient conditions (high COPs can only be achieved if 
a system’s thermalization time is short compared to the characteristic cycle 
time). 

The thermal conductivity of dielectric salts is on the order of ICO 2 W/(K 
■ m) at 0.1 K [13]. As a consequence, the thermal path in the salt must be 
minimized, and crystals of these materials are often grown around heat busses 
made of bundles of highly conductive metal wires. Figure 17 shows such a 
heat bus made from gold wires, which has been used for the construction of 
an FAA salt pill. Figure 18 shows a different design which has been devised 
at NASA GSFC for CPA salt pills. 

To optimize the design, the size of the metal structure providing the 
thermal conductivity has to be traded off against the lost salt volume. I 11 a 
simplified thermal model of a refrigerant pill, shown in Figure 19, the steady 
state heat flux through the heat bus into the refrigerant at temperature T r 
can be calculated from the following differential equation: 


d 2 T(x) 

dx 2 


a 


0 


(T(x)-T r ), 


(26) 


A • A hb 

where a characterizes the boundary resistance between the heat bus and the 
refrigerant, [3 is the contact area between the heat bus and the refrigerant 
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Figure 17: A heat bus made from a bundle of gold wires around which FAA crystals were 
grown for the salt pill of the ADR for the Micro-X sounding rocket telescope [85]. The 
gold wires ensure good thermal contact between the salt crystals and the two heat busses 
of the salt pill. 


per unit length, A the thermal conductivity of the heat bus material and Ahb 
the cross section of the heat bus. 

With the boundary conditions that T(x = 0) = T co i d and Q(x = l) = 0, 
one obtains the following solution for the temperature profile in the heat bus: 


where 


T{x) =T r + 


Tcoid — T, 

1 + e 


R 


2 ml 


(e K 


^ 2 ml ' g—mx 


). 


(27) 


m = 


a ■ P 


A • Ahb 

The total heat flux from the heat bus into the refrigerant is 


(28) 


HI 1 

Q(x = 0) = —A • Ahb • ( ^7 


X=0 


= -A- A 


HB 


m ■ ( T co id — Tr ) 


1 + e 


2 ml 


(1 — e 2mi ). (29) 


The total amount of heat Q that can be absorbed at the temperature T cold 
is proportional to the product of the amount of refrigerant (A tot — Ahb) • l 
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Figure 18: A copper heat bus for a CPA salt pill which was built at the NASA Goddard 
Space Flight Center. As CPA is not as corrosive as FAA, the heat bus for CPA salt pills 
can be constructed from copper. The copper part was manufactured with an Electrical 
Discharge Machining (EDM) process. A stainless steel sleeve is welded around the copper 
heat bus before the crystals are grown. 


and the total heat flux Q into the refrigerant for a given combination of T coW 
and T r . The maximum of Q can be obtained by differentiation with respect 

to Afffj. 

The fluorides and oxides which are used as magnetocaloric refrigerants 
at higher temperatures (such as GGG or GLF) can often be grown as large 
single crystals. The thermal conductivity is on the order of 10 W/(K • m) 
at 3 K and 100 W/(K • m) at 10 K [86], and no special measures need to be 
taken to enhance the thermal contact to the refrigerant. Very often, fluoride 
and oxide refrigerants are simply sealed in sintered or powdered form into 
hermetic containers with some exchange gas. 
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(a) ( b ) 

Figure 19: A simplified model of a refrigerant pill. In this model, the refrigerant has 
infinite thermal conductivity in the radial direction, and zero thermal conductivity in the 
axial direction. Axial heat conduction occurs only in the heat bus. In a real refrigerant 
pill, a large number of metal wires is used instead of one single rod to approximate this 
idealized model. 


5. CONCLUSION 

In cryogenics, magnetic refrigeration remains the most suitable technol- 
ogy for many applications. In most cases, magnetic refrigeration systems 
are optimized such that the volume of the refrigerant is minimal. In this 
paper, we have reviewed the most important properties of magnetic refriger- 
ants which influence cooling power density, followed by a description of the 
magnetocaloric materials most commonly used in low temperature applica- 
tions (paramagnetic salts, rare-earth oxides and fluorides, and Laves-phase 
materials) . 
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